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Observational tests for the homogeneity of the Universe on large scales are re- 
viewed. Assuming the Cosmological Principle we then estimate cosmological pa- 
rameters by joint analysis of the Cosmic Microwave Background, Supernovae la, 
peculiar velocities, cluster abundance and redshift surveys. Our results are con- 
sistent with results obtained by other groups, suggesting a best- fit FRW Universe 
with Q,ra = 1 — A « 0.3 and Ha « 75 km/sec/Mpc. We point out some potential 
problems with this currently popular model. 



1 Introduction 

The Cosmological Principle was first adopted when observational cosmol- 
ogy was in its infancy; it was then little more than a conjecture, embody- 
ing 'Occam's razor' for the simplest possible model. Observations could not 
then probe to significant redshifts, the 'dark matter' problem was not well- 
established and the Cosmic Microwave Background (CMB) and the X-Ray 
Background (XRB) were still unknown. If the Cosmological Principle turned 
out to be invalid then the consequences to our understanding of cosmology 
would be dramatic, for example the conventional way of interpreting the age 
of the Universe, its geometry and matter content would have to be revised. 
Therefore it is important to revisit this underlying assumption in the light of 
new galaxy surveys and measurements of the background radiations. 

Like with any other idea about the physical world, we cannot prove a 
model, but only falsify it. Proving the homogeneity of the Universe is in 
particular difficult as we observe the Universe from one point in space, and 
we can only deduce isotropy directly. The practical methodology we adopt is 
to assume homogeneity and to assess the level of fluctuations relative to the 
mean, and hence to test for consistency with the underlying hypothesis. If 
the assumption of homogeneity turns out to be wrong, then there are numer- 
ous possibilities for inhomogeneous models, and each of them must be tested 
against the observations. 

Here we examine the degree of smoothness with scale by considering red- 
shift and peculiar velocities surveys, radio-sources, the XRB, the Ly-a forest, 
and the CMB. We discuss some inhomogeneous models and show that a frac- 
tal model on large scales is highly improbable. Assuming an FRW metric 
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we evaluate the 'best fit Universe' by performing a joint analysis of cosmic 
probes. 

2 Cosmological Principle(s) 

Cosmological Principles were stated over different periods in human history 
based on philosophical and aesthetic considerations rather than on funda- 
mental physical laws. Rudnicki (1995) summarized some of these principles 
in modern-day language: 

• The Ancient Indian: The Universe is infinite in space and time and is 
infinitely heterogeneous. 

• The Ancient Greek: Our Earth is the natural centre of the Universe. 

• The Copernican CP: The Universe as observed from any planet looks 
much the same. 

• The Generalized CP: The Universe is (roughly) homogeneous and 
isotropic. 

• The Perfect CP: The Universe is (roughly) homogeneous in space and 
time, and is isotropic in space. 

• The Anthropic Principle: A human being, as he/she is, can exist only 
in the Universe as it is. 

We note that the Ancient Indian principle can be viewed as a 'fractal 
model'. The Perfect CP led to the steady state model, which although more 
symmetric than the CP, was rejected on observational grounds. The An- 
thropic Principle is becoming popular again, e.g. in 'explaining' a non-zero 
cosmological constant. Our goal here is to quantify 'roughly' in the definition 
of the generalized CP, and to assess if one may assume safely the Friedmann- 
Robertson- Walker (FRW) metric of space-time. 

3 Probes of Smoothness 

3.1 The CMB 

The CMB is the strongest evidence for homogeneity. Ehlers, Garen and Sachs 
(1968) showed that by combining the CMB isotropy with the Copernican 
principle one can deduce homogeneity. More formally the EGS theorem (based 
on Liouville theorem) states that "If the fmidamental observers in a dust 
spacetime see an isotropic radiation field, then the spacetime is locally FRW" . 
The COBE measurements of temperature fluctuations AT/T = 10~^ on scales 
of 10° give via the Sachs Wolfe effect {AT/T = ^A(f)/(?) and Poisson equation 



lisbon3: submitted to World Scientific on February 1, 2008 2 



rms density fluctuations of ^ ^ 10 ^ on 1000 h ^ Mpc (e.g. Wu, Lahav & 
Rees 1999; see Fig 3 here), i.e. the deviations from a smooth Universe are 
tiny. 

3.2 Galaxy Redshift Surveys 

Figure 1 shows the distribution of galaxies in the ORS and IRAS redshift 
surveys. It is apparent that the distribution is highly clumpy, with the Super- 
galactic Plane seen in full glory. However, deeper surveys such as LCRS and 
2dFGRS (sec below) show that the fluctuations decline as the length-scales 
increase. Peebles (1993) has shown that the angular correlation functions for 
the Lick and APM surveys scale with magnitude as expected in a universe 
which approaches homogeneity on large scales. 

Multifibre technology now allows us to measure redshifts of millions of 
galaxies. Two major surveys are underway. The US Sloan Digital Sky Survey 
(SDSS) will measure redshifts to about 1 million galaxies over a quarter of the 
sky. The Anglo- Australian 2 degree Field Galaxy Redshift Survey (2dFGRS) 
survey will measure redshifts for 250,000 galaxies selected from the APM 
catalogue. About 150,000 2dF redshifts have been measured so far (as of 
March 2001). The median redshift of both the SDSS and the 2dFGRS galaxy 
redshift surveys is ^ ^ 0.1. While they can provide interesting estimates of 
the fluctuations on intermediate scales (e.g. Peacock et al. 2001; sec Fig 2), 
the problems of biasing, evolution and ii'-correction, would limit the ability 
of SDSS and 2dF to 'prove' the Cosmological Principle, (cf. the analysis of 
the ESO slice by Scaramella et al 1998 and Joyce et al. 1999). 

3.3 Peculiar Velocities 

Peculiar velocities are powerful as they probe directly the mass distribution 
(e.g. Dekel et al. 1999). Unfortunately, as distance measurements increase 
with distance, the scales probed are smaller than the interesting scale of transi- 
tion to homogeneity. Conflicting results on both the amplitude and coherence 
of the flow suggest that peculiar velocities cannot yet set strong constraints 
on the amplitude of fluctuations on scales of hundreds of Mpc's. Perhaps 
the most promising method for the future is the kinematic Sunyaev-Zeldovich 
effect which allows one to measure the peculiar velocities of clusters out to 
high redshift. 

The agreement between the CMB dipole and the dipole anisotropy of 
relatively nearby galaxies argues in favour of large scale homogeneity. The 
IRAS dipole (Strauss et al 1992, Webster et al 1998, Schmoldt et al 1999) 
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Figure 1. The distribution of galaxies projeeted on the sky in the IRAS and ORS samples. 
This is an AitofT projection in Supcrgalactic coordinates, with L = 90°, B = 0° (close to the 
Virgo cluster) in the centre of the map. Galaxies within 2000 km/sec are shown as circled 
crosses; galaxies between 2000 and 4000 km/sec arc indicated as crosses, and dots mark 
the positions of more distant objects. Here we include only catalogued galaxies, which is 
why the Zone of Avoidance is so prominent in these two figures. (Plot by M. Strauss, from 
Lahav et al. 2000). 



shows an apparent convergence of the dipole, with misaUgnment angle of only 
15°. Schmoldt ct al. (1999) claim that 2/3 of the dipole arises from within a 
40 Mpc, but again it is difficult to 'prove' convergence from catalogues of 
finite depth. 
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Figure 2. The distribution of galaxies in part of the 2dFGRS, drawn from a total of 141,402 
galaxies. The slices are 4 dcg thick, towards the Northern Galactic Pole (left) and towards 
the Southern Galactic Pole (right). Not all 2dF fields within the slice have been observed 
at this stage: hence there arc weak variations of the density of sampling. The image reveals 
a wealth of structure, including superclusters and voids, but the similarity between the 
two slices suggests that on large scales the universe is isotropic and homogeneous, (from 
Peacock et al. 2001). 



3.4 Radio Sources 

Radio sources in surveys have typical median redshift ^ ~ 1, and hence are 
useful probes of clustering at high redshift. Unfortunately, it is difficult to 
obtain distance information from these surveys: the radio luminosity function 
is very broad, and it is difficult to measure optical redshifts of distant radio 
sources. Earlier studies claimed that the distribution of radio sources supports 
the 'Cosmological Principle'. However, the wide range in intrinsic luminosities 
of radio sources would dilute any clustering when projected on the sky. Recent 
analyses of new deep radio surveys (e.g. FIRST) suggest that radio sources 
are actually clustered at least as strongly as local optical galaxies (e.g. Cress 
et al. 1996; Magliocchetti et al. 1998). Nevertheless, on the very large scales 
the distribution of radio sources seems nearly isotropic. Comparison of the 
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measured quadrupole in a radio sample in the Green Bank and Parkes-MIT- 
NRAO 4.85 GHz surveys to the theoreticaUy predicted ones (Baleisis et ah 
1998) offers a crude estimate of the fluctuations on scales A ~ 600/i~^ Mpc. 
The derived amplitudes are shown in Figure 3 for the two assumed Cold 
Dark Matter (CDM) models. Given the problems of catalogue matching and 
shot-noise, these points should be interpreted at best as 'upper limits', not as 
detections. 

3.5 The XRB 

Although discovered in 1962, the origin of the X-ray Background (XRB) is 
still unknown, but is likely to be due to sources at high redshift (for review see 
Boldt 1987; Fabian & Barcons 1992). The XRB sources are probably located 
at redshift z < b, making them convenient tracers of the mass distribution 
on scales intermediate between those in the CMB as probed by COBE, and 
those probed by optical and IRAS redshift surveys (sec Figure 3). 

The interpretation of the results depends somewhat on the nature of the 
X-ray sources and their evolution. By comparing the predicted multipoles to 
those observed by HEAOl (Lahav et al. 1997; Treycr ct al. 1998; Scharf et 
al. 2000) we estimate the amplitude of fluctuations for an assumed shape of 
the density fluctuations (e.g. CDM models). Figure 3 shows the amplitude 
of fluctuations derived at the effective scale A ^ 60G/i^^ Mpc probed by 
the XRB. The observed fluctuations in the XRB are roughly as expected 
from interpolating between the local galaxy surveys and the COBE CMB 
experiment. The rms fluctuations ^ on a scale of 600/i~^Mpc are less than 
0.2 %. 

3.6 The Lyman-a Forest 

The Lyman-a forest reflects the neutral hydrogen distribution and therefore 
is likely to be a more direct trace of the mass distribution than galaxies are. 
Unlike galaxy surveys which are limited to the low redshift Universe, the for- 
est spans a large redshift interval, typically 1.8 < 2 < 4, corresponding to 
comoving interval of ~ 600 ft." ^ Mpc. Also, observations of the forest are not 
contaminated by complex selection effects such as those inherent in galaxy 
surveys. It has been suggested qualitatively by Davis (1997) that the absence 
of big voids in the distribution of Lyman-a absorbers is inconsistent with the 
fractal model. Furthermore, all lines-of-sight towards quasars look statisti- 
cally similar. Nusser & Lahav (2000) predicted the distribution of the flux in 
Lyman-a observations in a speciflc truncated fractal-like model. They found 
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Figure 3. A compilation of density fluctuations on different scales from various observa- 
tions: a galaxy survey, deep radio surveys, the X-ray Background and Cosmic Microwave 
Background experiments. The measurements arc compared with two popular Cold Dark 
Matter models (with normalization erg = I and shape parameters F = 0.2 and 0.5). The 
Figure shows mean-square density fluctuations (^)^ oc k^P{k), where k = 1/A is the 
wavenumber and P{k) is the power-spectrum of fluctuations. The open squares at small 
scales are estimates from the APM galaxy catalogue (Baugh & Efstathiou 1994). The elon- 
gated 'boxes' at large scales represent the COBE 4-yr (on the right) and Tenerife (on the 
left) CMB measurements (Gawiser & Silk 1998). The solid triangles and crosses represent 
amplitudes derived from the quadrupole of radio sources (Baleisis et al. 1998) and the 
quadrupole of the XRB (Lahav et al. 1997; Treyer et al. 1998). Each pair of estimates 
corresponds to assumed shape of the two CDM models. (A compilation from Wu, Lahav & 
Rees 1999). 



that indeed in this model there are too many voids compared with the ob- 
servations and conventional (CDM-like) models for structure formation. This 
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too supports the common view that on large scales the Universe is liornoge- 
neous. Another test for isotropy, based on the distribution of Supernovae la 
out to redshift 2; » 1 is described in Kolatt & Lahav (2001). 

4 Is the Universe a Fractal ? 

The question of whether the Universe is isotropic and homogeneous on large 
scales can also be phrased in terms of the fractal structure of the Universe. 
A fractal is a geometric shape that is not homogeneous, yet preserves the 
property that each part is a reduced-scale version of the whole. If the mat- 
ter in the Universe were actually distributed like a pure fractal on all scales 
then the Cosmological Principle would be invalid, and the standard model 
in trouble. As shown in Figure 3 current data already strongly constrain 
any non-uniformities in the galaxy distribution (as well as the overall mass 
distribution) on scales > 300/i~^Mpc. 

If we count, for each galaxy, the number of galaxies within a distance R 
from it, and call the average number obtained A'^(< R), then the distribution 
is said to be a fractal of correlation dimension D2 ii N{< R) oc i?^^. Of 
course D2 may be 3, in which case the distribution is homogeneous rather 
than fractal. In the pure fractal model this power law holds for all scales of 
R. 

The fractal proponents (Pietronero et al. 1997) have estimated D2 « 2 
for all scales up to ~ 500 Mpc, whereas other groups have obtained scale- 
dependent values (for review see Wu et al. 1999 and references therein). 

Estimates of D2 from the CMB and the XRB are consistent with D2 = 3 
to within 10"'^ on the very large scales (Peebles 1993; Wu et al. 1999). While 
we reject the pure fractal model in this review, the performance of CDM-like 
models of fluctuations on large scales have yet to be tested without assuming 
homogeneity a priori. On scales below, say, 30/i~^Mpc, the fractal nature 
of clustering implies that one has to exercise caution when using statistical 
methods which assume homogeneity (e.g. in deriving cosmological parame- 
ters). We emphasize that we only considered one 'alternative' here, which is 
the pure fractal model where D2 is a constant on all scales. 

5 More Realistic Inhomogeneous Models 

As the Universe appears clumpy on small scales it is clear that assuming the 
Cosmological Principle and the FRW metric is only an approximation, and 
one has to average carefully the density in Newtonian Cosmology (Buchert 

6 Ehlers 1997). Several models in which the matter in clumpy (e.g. 'Swiss 
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cheese' and voids) have been proposed (e.g. Zeldovich 1964; Krasinski 1997; 
Kantowski 1998; Dyer & Roeder 1973; Holz & Wald 1998; Celerier 1999; 
Tomita 1999). For example, if the hne-of-sight to a distant object is 'empty' it 
resuhs in a gravitational lensing de-magnification of the object. This modifies 
the FRW luminosity-distance relation, with a clumping factor as another free 
parameter. When applied to a sample of SNe la the density parameter of 
the Universe fim could be underestimated if FRW is used (Kantowski 1998; 
Perlmutter et al. 1999). Metcalf & Silk (1999) pointed out that this effect 
can be used as a test for the nature of the dark matter, i.e. to test if it is 
smooth or clumpy. 

6 Cosmological Parameters from a Joint Analysis: a Cosmic 
Heirmony ? 

A simultaneous analysis of the constraints placed on cosmological parameters 
by different kinds of data is essential because each probe (e.g. CMB, SNe 
la, redshift surveys, cluster abundance and peculiar velocities) typically con- 
strains a different combination of parameters. By performing joint likelihood 
analyses, one can overcome intrinsic degeneracies inherent in any single analy- 
sis and so estimate fundamental parameters much more accurately. The com- 
parison of constraints can also provide a test for the validity of the assumed 
cosmological model or, alternatively, a revised evaluation of the systematic 
errors in one or all of the data sets. Recent papers that combine informa- 
tion from several data sets simultaneously include Webster et al. (1998); 
Lineweavcr (1998); Gawiser & Silk (1998), Bridle et al. (1999, 2001), Eisen- 
stein, Hu & Tegmark 1999; Efstathiou et al. 1999; and Bahcall et al. (1999). 

While joint Likelihood analyses employing both CMB and LSS data are 
allowing more accurate estimates of cosmological parameters, they involve 
various subtle statistical issues: 

• There is the uncertainty that a sample does not represent a typical patch 
of the FRW Universe to yield reliable global cosmological parameters. 

• The choice of the model parameter space is somewhat arbitrary. 

• One commonly solves for the probability for the data given a model (e.g. 
using a Likelihood function), while in the Bayesian framework this should 
be modified by the prior for the model and its parameters. 

• If one is interested in a small set of parameters, should one marginalise 
over all the remaining parameters, rather than fix them at certain (some- 
what ad-hoc) values ? 
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• The 'topology' of the Likelihood contours may not be simple. It is helpful 
when the Likelihood contours of different probes 'cross' each other to yield 
a global maximum (e.g. in the case of CMB and SNe), but in other cases 
they may yield distinct separate 'mountains', and the joint maximum 
Likelihood may lie in a 'valley'. 

• Different probes might be spatially correlated, i.e. not necessarily inde- 
pendent. 

• What weight should one give to each data set ? 

In a long term collaboration in Cambridge (Bridle et al. 1999, 2001; 
Efstathiou et al. 1999; Lahav et al. 2000) we have compared and combined in 
a self-consistent way the most powerful cosmic probes: CMB, galaxy redshift 
surveys, galaxy cluster number counts, type la Supernovae and galaxy peculiar 
velocities. These analyses illustrate the power of combining different data sets 
for constraining the fundamental parameters of the Universe. Our analysis 
suggests, in agreement with studies by other groups, that we live in a flat 
accelerating Universe, with comparable amounts of dark matter and 'vacuum 
energy' (cosmological constant). We have also addressed recently (Lahav et 
al. 2000; Lahav 2001) the issue of combining diS^erent data sets, which may 
suffer different systematic and random errors. We generalised the standard 
procedure of combining likelihood functions by allowing freedom in the relative 
weights of various probes. This is done by including in the joint likelihood 
function a set of 'Hyper-Parameters', which are dealt with using Bayesian 
considerations. The resulting algorithm, which assumes imiform priors on 
the logarithm of the Hyper-Parameters, is simple to implement. Here we 
show some examples of and results from the joint analysis. First combining 
two CMB data sets, and then combining CMB, Supernovae la and peculiar 
velocities. 

7 Combining the CMB Boomerang and Maxima Data 

The recent Boomerang (hereafter B; de Bernardis et al. 2000) and Max- 
ima (hereafter M; Hanany et al. 2000) CMB anisotropy measurements 
yielded high-quality angular power spectra Ci over the spherical harmonics 
400 / ^ 800. An important factor in interpreting the data is the calibra- 
tion error. The experimental papers quote calibration errors of 10% and 4% 
(1-sigma in AT/T) for B and M, respectively. The measurements (with B 
data corrected upward by 10%, and M data corrected downward by 4 %) 
are shown in Figure 4, and they indicate a well defined first acoustic peak 
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at I ~ 200, with less convincing second and third peaks at higher harmon- 
ics. These measurements favour (under certain assumptions) a flat universe, 
spectral index n = 1 and baryon density fl^h^ ~ 0.03 (e.g. Jaffe et al. 2000), 
which is about 2-sigma higher than the Big-Bang Nucleosynthesis (BBN) value 
n^h^ ~ 0.0190 ± 0.0018 (95 % CL; Buries et al. 2000). Note that the recent 
CBI result (Padin et al. 2000) gives a higher power (at I ~ 600) relative to 
B&M. Jaffe et al. (2000) fitted models after combining the B& M data sets 
into one set. We have taken a different approach for joint analysis of these 
two data sets by utilising the 'Hyper-Parameters' (Lahav 2001). 

For simplicity, we restrict ourselves to a very limited set of cosmological 
models. We obtain theoretical CMB power-spectra using the CMBFAST and 
CAME codes (Slejak & Zaldarriaga 1996; Lewis, Challinor & Lasenby 2000). 
We assume that CMB fluctuations arise from adiabatic initial conditions with 
Cold Dark Matter (CDM) and negligible tensor component, in a flat Universe 
with fim = 0.3, A = A/(3iJ^) = 1 - fi^n = 0.7, n = 1, Qrms = IS/itK and 
r^b^^ = 0.03. This choice is motivated by numerous other studies which 
combined CMB data with other cosmological probes (e.g. Jaffe et al. 2000, 
Bridle et al. 2000; Hu et al. 2000; see also the next section). We then 
investigate the constraints on the remaining parameter, the dimensionless 
Hubble constant, h = ffo/(100 km/sec/Mpc). Increasing h decreases the 
height of the first acoustic peak, and makes few other significant changes to 
the angular power spectrum (e.g. Hu et al. 2000). The range in h investigated 
here is (0.5 < h < 1.1). 

The calibration of the data brings the two data sets to much better agree- 
ment. In fact, in this case the standard joint ^J^d flic Hyper Parameters 
give the same result, h = 0.79, with slightly smaller error bars in the HPs 
case (±0.04; 95% CL). We also tried the BBN value fl^h^ = 0.019 which gives 
gives much poorer than the value flbh"^ = 0.03 (e.g. Jaffe et al. 2000 and 
others). 

The best fit Hubble constant is i?o = 79 ± 4 km/sec/Mpc (95% CL, 

random errors only) for a fixed flat CDM = 1 — A = 0.3 model with 
n = 1, Qrms = 18/xK and Q\)h^ = 0.03. We note that if more cosmological 
parameters are left free and then marginalised over, the error in h would 
typically be larger. 

This combination of fim and Hq corresponds gives for the age of the Uni- 
verse 11.9 Gyr. Our derived Hq is slightly higher but still consistent with the 
'flnal result' of Ha from Cepheids and other distance indicators (Freedman et 
al. 2000) Hq = 72± (3)r-± (7)s km/sec/Mpc (1-sigma random and systematic 
errors). 



lisbon3: submitted to World Scientific on February 1, 2008 11 



Figure 4. The Boomerang -^^ data against spherical harmonic I (top panel, calibrated by 
1.10) and Maxima data (bottom panel, calibrated by 0.96). The line in each panel is for a 
A-CDM model with n = l,Qm = l- A = 0.3, Qb/*^ = 0.03, Qrms = 18nK, and our 'best 
fit' h = 0.8. 



8 Combining CMB, Supernovae la and Peculiar Velocities 

A recent study (Bridle et al. 2001) is an example of combining 3 different data 
sets. We compared and combined likelihood functions for the matter density 
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parameter fiiji, the Hubble constant h, and the normahzation (in tc^rnis of 
the variance in the mass density field measured in an 8h~^ Mpc radius sphere) 
from peculiar velocities, CMB (including the Boomerang and Maxima data) 
and type la Supernovae. These three data sets directly probe the mass in the 
Universe, without the need to relate the galaxy distribution to the underlying 
mass via a "biasing" relation. 

Our analysis assumes a flat ACDM cosmology with a scale-invariant adi- 
abatic initial power spectrum and baryonic fraction as inferred from big-bang 
nucleosynthesis. We find that all three data sets agree well, overlapping sig- 
nificantly at the 2(7 level. This therefore justifies a joint analysis, in which 
we find a best fit model and 95% confidence limits of = 0.28(0.17,0.39), 
h = 0.74(0.64,0.86), and as = 1.17(0.98,1.37). In terms of the natural 
parameter combinations for these data agQ^'^'^ = 0.54(0.40,0.73), ilmh = 
0.21 (0.16,0.27). Also for the best fit point, Qrms-ps = 19.7/xK and the age of 
the Universe is 13.0 Gyr. Figure 5 illustrates these results. 

This is quite in agreement from results form cluster abundance Om'^'^crg = 
0.5 ± 0.1. (Eke et al. 1998). By combining the abundance of clusters with 
the CMB and IRAS Bridle et al. (1999) found Q^ = l-X = 0.36, h = 0.54, 
(Tg ~ 0.74, and IRAS biasing parameter biras = 1-08 (with error bars similar 
to those above). 

9 Discussion 

Analysis of the CMB, the XRB, radio sources and the Lyman-a which probe 
scales of ~ 100 — 1000 ft. ~^ Mpc strongly support the Cosmological Principle 
of homogeneity and isotropy. They rule out a pure fractal model. However, 
there is a need for more realistic inhomogeneous models for the small scales. 
This is in particular important for understanding the validity of cosmological 
parameters obtained within the standard FRW cosmology. 

Joint analyses of the CMB, IRAS, SNe, cluster abundance and peculiar 
velocities suggests = 1 — A « 0.3. The measurement of the Hubble constant 
from Cepheids and from the CMB suggests Hq ps 75 km/sec/Mpc. While 
this is now a popular model there are potential problems with this set of 
parameters. There is no simple theoretical explanation why the present epoch 
contributions to matter Om and 'dark energy' (A) are nearly equal. Also, for 
the above model the age on the Universe is ~ 13 Gyr which is uncomfortably 
close to some estimates for the age of the Globular Clusters. 

These issues will not doubt be revisited soon with larger and more accu- 
rate data sets. We will soon be able to map the fluctuations with scale and 
epoch, and to analyze jointly redshift surveys (2dF, SDSS) and CMB (MAP, 
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Figure 5. The 1-dimensional marginalised likelihood distributions from the joint PV, 
CMB and SN likelihood function. The main results are shown by the solid lines, which 
use the whole CMB data compilation, PV (marginalised over the field velocity dispersion 
(Tv) and SN. The dotted lines show the likelihood functions when PV, SNe and just the 
pre-BM CMB data are used. The dashed line is the result when the CMB data is just 
COBE + BOOMERANG + MAXIMA-1. The result (using all the CMB data) when the 
uncorrelated velocity dispersion term is not included in the PV analysis (ctv = 0) is shown 
by the dot-dashed line. (From Bridle etal. 2001.) 



Planck) data. These high quahty data sets will allow us to study a wider 
range of models and parameters. 
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